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We continue the recent study of the possibility of constructing a consistent infrared modifi- 
cation of gravity by treating the vierbein and connection as independent dynamical fields. 
We present the generalized Fierz-Pauli equation that governs the propagation of a massive 
spin-2 mode in a model of this sort in the backgrounds of arbitrary torsionless Einstein 
manifolds. We show explicitly that the number of propagating degrees of freedom in these 
backgrounds remains the same as in flat space-time. This generalizes the recent result that 
the Boulware-Deser phenomenon does not occur in de Sitter and anti-de Sitter backgrounds. 
We find that, at least for weakly curved backgrounds, there are no ghosts in the model. We 
also briefly discuss the interaction of sources in flat background. 
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1 Introduction and summary 



The possibility that gravity may get modified at large distances attracts considerable in- 
terest, which is motivated, in particular, by the accelerated expansion of the Universe. 
Deforming General Relativity in the infrared is, however, not at all easy. The problems 
one has to worry about are well understood in the context of the Fierz-Pauli theory pQ of 
massive graviton. Already at the linearized level about Minkowski background, the Fierz- 
Pauli gravity exhibits the van Dam-Veitman-Zakharov discontinuity [21 [3], which, taken at 
face value, implies that the theory deviates from General Relativity even at short distances. 
Below the Vainshtein radius, however, non-linear effects become important [3j and, in fact, 
these effects may cure the van Dam-Veltman-Zakharov problem [5]. The disaster occurs in 
curved backgrounds, where an extra, Boulware-Deser propagating mode shows up [HO [8], 
over and beyond the five modes of massive graviton. This extra mode has ghost kinetic 
term and makes the theory unacceptable. 

There are various theories that pretend to be free of at least some of the above prob- 
lems. These include theories which deform General Relativity in backgrounds other than 
Minkowski [91 [TOl lllj . theories with extra dimensions [12[ [T3"l HM \TE[ [T6] (for a review 
see Ref. |17| ) and theories with broken Lorentz-invariance |18[ [191 [20], [2Tj [22] (for a re- 
view see Ref. [23]). Yet another option is to consider theories whose independent fields 
are both vierbein and connection and whose Lagrangians contain terms quadratic in cur- 
vature as well as mass terms for the torsion field. Some theories of the latter type have 
both massless and massive spin-2 modes in the spectrum about Minkowski background 
and nevertheless are free of pathologies (ghosts or tachyons) at the linearized level in this 
background [25 E3 EBl EZ! • 

It is worth noting that theories with dynamical vierbein and connection may be viewed as 
gauge theories % with spontaneously broken local Lorentz invariance [301131] . This viewpoint 
opens up a possibility of unification of gravity with the Yang-Mills theory [321 [331 [M]. We 
will not pursue this approach and concentrate on degrees of freedom describing dynamical 
space-time geometry itself. 

Once a geometrical theory has a massive spin-2 mode at the linearized level in Minkowski 
background, a question arises of whether or not this theory has the Boulware-Deser mode 
in curved backgrounds. Recently, this question has been addressed |35] in the context of 

2 Vierbein and connection as gauge fields of the Poincare group have been introduced by Kibble [2S] and 
studied extensively by Hehl and collaborators. For a historical review and references to earlier work see 
Ref. |29| . Our point here is to study the possibility of a consistent infrared modification of gravity. 



2 



one of the models of Refs. \24\ [25l [26l I27j . The analysis in Ref. |35] has been restricted 
to de Sitter and anti-de Sitter backgrounds. The result has been encouraging, as it has 
been found that there are no Boulware-Deser modes in these backgrounds. In this paper 
we extend the analysis to arbitrary Einstein backgrounds with vanishing torsion and show 
that the Boulware-Deser mode does not appear in this more general case as well. If the 
curvature of the background is sufficiently small, propagating modes are neither ghosts nor 
tachyons. Hence, our analysis supports the conjecture of Ref. [35] that the model is healthy 
at least for sufficiently weak fields. 

As a by-product, we identify the field that reduces to the massive spin-2 field as the 
curvature of the background is switched off, and obtain the equation for this field that serves 
as a generalization of the Fierz-Pauli equation to arbitrary Einstein backgrounds. 

In the end of this paper we consider the interation of sources in our model at the 
linearized level in Minkowski background. We rederive the result of Refs. [25] that within 
the appropriate range of parameters, there are no ghosts or tachyons. The interaction 
between symmetric and conserved energy-momentum tensors is mediated by both massless 
and massive spin-2 fields, the relative strength depending on the parameters of the model. 
In this sense our model is indeed an infrared-modified gravity. The interaction due to 
massive spin-2 field does exhibit the van Dam-Veltman-Zakharov discontinuity. We leave 
for future study the questions of whether the Vainshtein phenomenon occurs in our model, 
and if so, whether it cures the van Dam-Veltman-Zakharov problem. 

The paper is organized as follows. We introduce the model, its field equations and 
the Einstein manifold solutions in section [2j In section [3] we discuss the behaviour of 
vector components of the torsion field in the Einstein backgrounds. The results of section [3] 
have been already obtained in Ref. [35]; we present them for completeness and later use. 
Section [4] contains the main results of this paper. In section [4~T1 we derive the equation that 
the massive tensor field perturbations obey in the Einstein backgrounds. This equation 
may be viewed as a generalization of the Fierz-Pauli equation. We then proceed to show 
that in arbitrary Einstein backgrounds, this equation describes five propagating modes, the 
right number for massive spin-2 field. We begin in section 14.21 by counting the number 
of constraints — equations that contain time derivatives of at most first order — among 
the ten equations constituting the generalized Fierz-Pauli system. We find that there are 
five such constraints, which suggests that there are indeed five propagating modes. In 
section 14.31 we make use of the Stiickelberg formalism to confirm this result and show that 
at least in weakly curved backgrounds, neither of the propagating modes is a ghost. We 
note, however, that the ghost problem may reappear in backgrounds of sufficiently high 



3 



curvature. Finally, in section [5] we study the linearized theory with sources in Minkowski 
background and obtain both the fields induced by the sources and the action that describes 
the interaction of the sources. Some of the results of this section are contained in Ref. [25J; 
however, our emphasis will be different. We conclude in section [6] by discussing directions 
for future studies. 

2 The model 
2.1 Action 

We follow the notations of Refs. [2i\ |2"51 1251 E5] and denote the vierbein by and the 
connection by Aij^ = —Aji^, where \i = 0,1,2,3 and i,j = 0,1,2,3 are the space-time 
and tangent space indices, respectively. We often use the tangent space basis, in which the 
indices are raised and lowered by the Minkowski metric rjij. The connection can be viewed 
as an 0(1,3) gauge field. It is conveniently decomposed as follows, 

A-ijk = -^ijfJt^f- — 2 i^ijk T^jik 2~fcr/' ~i~ Cijfc Cjik Cfcij) ; (1) 

where 

is constructed from vierbein and = —Tn-j is the torsion tensor. The latter can in turn 
be decomposed into its irreducible components under the local 0(1, 3) group, 

T ijh = -(t ijh - t ikj ) + -( W - VikVj ) + s ijkl a l (2) 

where the field is symmetric with respect to the interchange of i and j and satisfies the 
cyclic and trace identities, 

tijk ~\~ tjki ~i~ tkij — 0, V^^ijk — 0> V ^ijk — 

The 24 independent components of break up into 4 components of Vi, 4 components of 
di and 16 independent components of 

The curvature, as usual in gauge theories, is defined by 

The model studied in this paper, as well as in Ref. [35], is defined by the action 

S = [ (fxe (L F + L T ) , 
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where e = det e^, 



and 



with 



L F = C1 F + c 2 + csFijFV + c^FijF^ + c 5 F 2 + c 6 {e ijkl F^ 



L 1: - a [ t ijk t ijk - vy + ^did' 



kl\2 



F jl = r, ik F i3kU F = r ] i k F jk , e ■ F = e ljkl F^ kl . (3) 

Here c\,...,cq are "coupling constants" obeying, apart from sign restrictions, the only 
condition 

c 3 + c 4 + 3c 5 = . 

In what follows, three combinations of these parameters will be used, 

2 

a = a + -ci (4) 
A = -±. (5) 

DCl 

x = 2A + — (6) 

2c 5 

By appearance, the term Lp has the form of kinetic term for the connection (plus the 
cosmological constant term C2), while Lt is torsion mass term. 

We note in passing that the Lagrangian Lj- does not explicitly break local 0(1,3) in- 
variance, so the entire action is invariant under both local frame rotations and general 
coordinate transformations. For C2 = 0, the model admits Minkowski space-time as a solu- 
tion of the field equations. In that case, the local 0(1, 3) invariance is spontaneously broken 
by the background value of the vierbein field, cf. Refs. [301 EH E21 1331 134"] . 

The model with C2 = is free of ghosts and tachyons in Minkowski background provided 
the parameters satisfy certain inequalities [241 125} [26 | I27j. which in our notations read 

ci > , c 5 < , c 6 > , a < , a > . (7) 

Non-vanishing value of C2 enables one to have de Sitter or anti-de Sitter solution with 
vanishing torsion in this model, with cosmological constant equal to A. In the latter case, 
the requirement of the absence of tachyons imposes one more condition [35], C5X > 0, i.e., 

a > -4Ac 5 . (8) 

Since C5 < 0, the latter condition is non-trivial for positive A. Once the above conditions 
are satisfied, the theory is healthy in de Sitter and anti-de Sitter backgrounds [35] . 
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2.2 Field equations 



There are two sets of field equations in our model. One consists of the gravitational field 
equations obtained by varying the action with respect to vierbein, 

ClFji + C$(F iF m j Fj j-fmn) C<l(F iFj m Fj j-^nm) ~i~ ^C^FjiF 

+2c 6 e kmnj F kmn i (e rpqs F r ™ s ) + (D k + v k )F ijk + H tj - \mA L F + L T ) = (9) 



where 



ijk 



{Ujk - tikj) - {Vij v k - VikVj) - -^ijkia 1 



tt rp rpmn ^LT 1 TP mn 



and Di is the covariant derivative with respect to the connection Aj_j k . Note that these 
equations have both symmetric and antisymmetric parts. 

By varying the action with respect to the connection A.^^ one finds another set of 
equations, 

c 3 \rf\D m + \v m )F^ m - rf k (D m + \v m )F im - (D l + + (LP + \v^)F lk ^ 

+c 4 U k (D m + \v m )F m i - rf k (D m + \v m )F m - (D* + -v^ + + ^)F ki \ 
+c 5 {^(D* + ^)F - rf k (D l + ^)f\ + 4c 6 {e^ m (An + |v m )(e • F) 

- (i fc [m n] + £ fc m n P a p ) { c 3 (r ? m F^ - if m F in ) + c 4 ( V im F^ - ^ F m ) 
+2c 5 7] im r}i n F + 2c 6 ^' fcm (e ■ F)| + H^ k = (10) 
where 

H ijk = -oi{t ki j - t kji ) + a(r] ki Vj - rj kj Vi) - —e ijk ia l 

Equations ([§]) and (fTUj) are not completely independent because of the Bianchi identity. In 
a theory with torsion, this identity reads 

D k Fiji m + T n k iF ijmn + cyclic (him) = (11) 

Contracting this identity one obtains 

- iDjF = 7*^1* i + (12) 

We will see that these identities provide useful constraints in the linearized theory. 
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2.3 Einstein backgrounds 

Let us consider torsion-free backgrounds. For vanishing torsion, the curvature tensor and 
its contractions reduce to the Riemann tensor, the Ricci tensor and the Ricci scalar, i^j-jy = 
Rijki, F{j = Rij and F = R, respectively. By inspecting eq. ([9]) one finds that it is satisfied 
for the Einstein manifolds. The Riemann tensor for these manifolds has the following form 

Rijki = MVikVjl ~ VUVjk) + Wijki, 

where the Weyl tensor Wijki has all symmetries of the Riemann tensor and is traceless in 
all pairs of indices. One then has 

Rij = SArjij , R = 12A . 

Using these properties one finds that eq. © reduces to the relationship © between A and 
coupling constants. Equation (jlOh is satisfied for the Einstein manifolds automatically. 

An important property of the Weyl tensor of the Einstein manifolds follows from the 
Bianchi identity, namely 

VW ti H = 

We will repeatedly make use of this property in what follows. 

3 Linearized theory in Einstein backgrounds: pseudovector 
di and vector Vi 

One of the main purposes of this paper is to study field perturbations about general torsion- 
free Einstein backgrounds. The analysis of pseudovector aj and vector V{ has been performed 
in Ref. [35] where is has been shown that the vector field does not have its own propagating 
modes, while the pseudovector field a% is a gradient, and its longitudinal part obeys the 
massive Klein-Gordon equation. These properties are exactly the same as in the theory 
about Minkowski background. For the sake of completeness and presentation of useful 
formulas, we recapitulate the analysis here. 

We do not use special notation for the background objects, unless there is a risk of an 
ambiguity; the subscript (1) refers to linearized perturbations. The torsion components 
vanish for our backgrounds, and we do not label their perturbations by the subscript (1). 

To study the fields Vi and dj, it suffices to consider certain combinations of the full 
equations Q and (jlOp . We begin with the antisymmetric part of the gravitational equation 
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([9]) whose complete form is 



+{D k + v k )F [ij]k + H m = . (13) 

By linearizing this equation about the Einstein background, one obtains 

(d - 4Ac 3 )F (1)bi] - V fe F (1)b , ]fc - (c 3 - c 4 )^[ m "lF (1)[mn] = . (14) 

where V denotes the covariant derivative with respect to the background metric. By lin- 
earizing Fij defined in (J3j) one finds the following explicit expression for the antisymmetric 
part 

F «w = ~L vkp ^ k (i5) 

Hence, eq. ([13]) is an algebraic equation for F^u^ whose only solution is 

F { im = (i6) 

This result simplifies all other equations. 

Now, let us write the curl of the torsion equation (|1U|) . Namely, we contract eq. (|1U|) 
with £ijki to find the complete curl equation, 

(C3 - c^D'F* - 12c„A(e ■ F) - 2 -e, m t*W + c lF «>) - S m (e ■ F) 

2 2 

-o( c 3 - c 4 )e ljkl v l F^ k - 2(c 3 F i; + c 4 i^V + -aa, = (17) 
o y 

Substituting (|16|) into eq. (|17|) and using the fact that e • F = 6V*Oj we obtain the following 
linearized equation for a/, 

8c 6 V,(V • a) - ^2Ac 5 + ~\ a, = (18) 

This equation shows that the pseudovector field is a gradient, ai = Via. Its longitudinal 
part obeys the Klein-Gordon equation 

where x is defined in ([6]) . The spin zero field a has healthy kinetic term and is not a tachyon 
provided that the inequality §8§ is satisfied. The mass of this field coincides with the flat 
space result when A = 0. 
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It is worth noting that using the explicit form of the right hand side of eq. (|15p and the 
fact that ai is a gradient, one obtains from eq. (|16p the following constraint 



V £fc[ mn ] = V[ m w n ] (20) 

This constraint already suggests that at least some components of the field Vi are not 
independent. 

In fact, the entire field vi is not dynamical by itself, as it can be expressed through the 
tensor tij k . To see this, one makes use of the trace of the torsion equation. For obtaining 
its general form, one contracts eq. (|10p with rf k and finds 

-3c 5 (DjFM - \d 1 F^ + (c 3 - c A )DjF^ - 2c 5 (vjF^ - \^f \ + (c 3 - c A )DjF^ 



+ l(c 3 - ct) Vi FM - 3c 5 f^F (jn) + l(c 3 - c,)f^F u 
1 



-(c 3 - c 4 )e ijnl aiF [jn] + Q^a^e ■ F) + -av 1 = (21) 



Now, one substitutes Fu^ = in the linearized version of this equation and finds 

DjF ij - -D'F) = xv l (22) 

The left hand side of eq. (|22p vanishes for vanishing torsion, so it is proportional to torsion 
at the linearized level. To obtain the explicit expression one makes use of the Bianchi 
identity (|12|) . Linearizing this identity in the Einstein background, one finds 

d% - ^Djf) = 2A Vj + \r kl w kl 3l 

1 J (1) 1 
so that eq. ([22]) becomes 

Vi = ^W ijH m (23) 

Hence, the vector field Vi is not an independently propagating field, exactly as in the flat 
space. 



4 Spin-2 mode in Einstein backgrounds 
4.1 Generalized Fierz-Pauli equation 

Let us now derive the equations for propagating tensor perturbations. To this end, we 
make use of the results of section [3] and write the liearized gravitational equation ^ in the 
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following form, 



c i " + ^ kp (i)ijk + ZcsWjmniF™ = (24) 



Hereafter we treat the field F(i)ij as well as the combination V k F^j_j k as symmetric with 
respect to the interchange of indices i and j, see eqs. (fT5j) and ([To]) . The expression for 
V fc -F(i)ijfe is, explicitly, 



2 

The remaining equation is the torsion equation (jlOp whose linearized version is 

1 



(25) 



^ (^2A + ^ [(^v* - W ) + 4t k[ij] ] = (26) 



Here Vi should be expressed in terms of tij k according to (|23|) . 

Equation (f26l) together with eq. (f22l) may be used to express i^-M in terms of 

(3V<F (1)ifc - 3V,F (1)ifc ) - (vikV l F m - VjkV l F {1)il/ 
+ {rHkV j F {1) - Vjk ViF {1) )] , (27) 



1 

4x 



where x is defined in (0). It is straightforward to check that eq. (f2"0"j) . with t^uu and -Uj 
expressed through i^nw, is identically satisfied. 



Equation (|27j) determines, in fact, the full tensor i^^. in terms of jPnvj. Indeed, due to 
the identities ([3]), one has 

= g (*i[jfc] + tj[ik\) (28) 



Substituting this into eq. ([25]) one obtains 



V Fj k i — — — 



4x 



6V 2 F (1)jfc - 3 (VjVF il)kl + W,F (1)H J - 3 (V k V l F (1)jl + VV fc F (1)j7 
+6 Vjk V m V n F (1)mn - 4r/ ifc V 2 F (1) + 4V i V fc F (1) ] (29) 

Thus, equation (124"]) becomes the equation for the field Fm^j, 

V k F ijk + ci yF^ i:j - ^T]ij F (i)J + 3c 5 W jmni F( 1)mn = , (30) 
where the first term is given by the right hand side of (|29p . 
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Equation ([30]) is a closed equation for the field -F(i)rp while the fields tijk and v,- t are 
expressed through Fn^j according to eqs. (p7|) . (f28|) and ([23]) . One can check that all 
equations of the linearized theory are satisfied provided that the field i^nvj obeys eq. (f3"Uj) . 



Equation (|30p does not look similar to the Fierz-Pauli equation yet. To write it in a 
more familiar way, let us introduce the field 

u ij = F (i)ij ~ QVijF(l) 

Then eq. ([30]) becomes 

V 2 Uij - V k ViU kj -V k VjU ki + ViVjU + raj (v k V l u H - V 2 u) + 6A ( Uij - -r^u 



- {2K + ^\ (u^ - THju) + (l - ^) W Uk] u lk = (31) 



This equation reduces to the Fierz-Pauli equation in Minkowski background, with the mass 
of the spin-2 field given by 

™ 2 = £ ■ < 32 » 

So, eq. (|3ip may be viewed as the generalization of the Fierz-Pauli equation to the Einstein 
backgrounds. 

Equation (|31|) has particularly simple form in de Sitter and anti-de Sitter backgrounds 
for which Wukj = 0. In that case the trace and divergence of this equation together with 
the equation itself imply that u%j has to be traceless and divergence free. We then obtain 
the Klein-Gordon equation with the mass given by 

M 2 = 4Afl + ^U " Cl 



3a / 3ac5 

in accord with Ref. [35]. This mass obeys the Higuchi bound [36] provided the inequality 
jHl) is satisfied. 



A remark is in order. Besides the tensor mode propagating according to eq. (I3ip . there 
is of course a massless tensor mode. The latter corresponds to perturbations of the vierbein 
field with vanishing torsion, and propagates according to the linearized Einstein equations 
(written relative to an orthonormal basis of the background geometry), Rn)ij = 0. For this 
mode = 0, so eq. (f3T|) is trivially satisfied. 

4.2 Counting the number of constraints 

Let us show that out of the ten equations satisfied by ten components of Uij, five are con- 
straints which involve at most first order time derivatives. The remaining five involve at 
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most second order time derivatives. This suggests that the field Ujj describes five propagat- 
ing modes, the right number for massive spin-2 field. 

Obtaining four out of five constraints is straightforward. Indeed, the divergence of 
eq. (|3Tj) gives first-order equations, 

M 2 V\u l3 - mj u) = (l- ^) W jikl V k u a (33) 

where 

M 2 = M 2 - 2A 

These are obviously the four constraints. 

To see that there are actually five constraints, let us choose the coordinates in the 
background manifold such that the background metric components are go a = = g 0a , 
where a = 1, 2, 3 and g 00 = g™. Working in coordinate basis rather than orthonormal one 
which we have used so far, we perform the (3 + l)-decomposition of the field components 
u^ u as well as the field equations. It is straightforward to see that (00)- as well as (a0)- 
components of eq. (|3ip involve at most first order time derivatives. These are the four 
constraints corresponding to eq. f|33|) . The (o6)-components of eq. (|3ip — six equations - 
are superficially second order in time. However, one combination of the latter is in fact a 
constraint. To see this, let us write the trace of eq. (f3Tj) . 



u = 2xc 5 ?L WimknV i V n u mk (34) 

M A xc\ 

To show that no time derivatives higher than first order enter this equation, we write the 
right hand side of this equation in an expanded form, 

u = ^Z£L W = 2 -^^W 0M V 2 u ab + .... (35) 

where dots denote the terms which have at most one time derivative. Now, we substitute 
VgU a fc from eq. (f3Tj) into eq. ([35]) . One can show that the expression for V\u a b does not 
involve any term with more than one Vo- Thus, the substitution shows that eq. (135p is 
indeed another constraint reducing the number of propagating modes from ten to five. 



4.3 Stiickelberg treatment 



A simple way to isolate the dangerous degrees of freedom is to make use of the Stiickelberg 
trick. As an example, this trick enables one to see in rather straightforward manner how the 
van Dam-Veltman-Zakharov phenomenon and related effects emerge |37j in the Fierz-Pauli 
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theory in Minkowski background, and how the Boulware-Deser ghost mode appears [8] 
in that theory in curved backgrounds. Let us make use of the Stiickelberg trick to see that 
no Boulware-Deser mode is present in the theory with the field equation ([31 j) . 



The quadratic action which corresponds to eq. (|3Tj) is, up to an overall factor, 

S = Si nv + S m + Sy/ 



where 



<l l .r^\ -^V fc uyVfcU y + V k u ki V lU H - V k u lk V l u + ^V^/V"// 



-A (uiju* + \u 2 ^ - W Ukj u lk uV | (36) 
s m = — 7T / (ftx^g {uijU lj - u 2 ) (37) 



Here 



2 

S w = s ( d^xV^g W iklj u kl u ij (38) 



2xc 5 
s = 1 

a 

The part (|36p of the action is invariant under the gauge transformation 

Uij = Uij + ViQ + VjCi (39) 

while the terms (I37p and (|38h are not invariant. To implement the Stiickelberg procedure, 
one introduces new fields Uy, £j and and writes 

Altogether, there are now 15 fields. The theory is invariant under two gauge symmetries, 

Uij > Uij + V iQj + Vj£j , £j > £j C« 

and 

6 - 6 + , - - 2^ (40) 

Provided that all fields have kinetic terms quadratic in derivatives, these symmetries elim- 
inate 10 degrees of freedom, and there remain 5 propagating modes. It is not at all guar- 
anteed, however, that the change of variables (|39p does not introduce higher-derivative 
terms in the action. A counterexample is given by the Fierz-Pauli gravity in curved back- 
grounds: the kinetic term for the field (j) has four derivatives (there are terms like (n^) 2 ), 
and the sixths propagating mode, a ghost, appears in the spectrum [23]; this is the 
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way the Boulware-Deser mode is seen in the Stiickelberg formalism. Let us show that this 
phenomenon does not occur in our model. 

The action Si nv = Si nv (u) contains only the field u^- and is obviously second order in 
derivatives. The change of variables (|39p gives rise to the derivative terms in the mass part 
of the action, 

M 2 f 

S m (u, £, 0) = -— J d^x^-g {(Vrf; - V^)(V^ - V j C) ~ 2uV 2 cf> - 3AV^V> + . . . } 

(41) 

where the contribution proportional to A appears in the process of integration by parts, and 
dots denote terms with less than two derivatives. The term (|38p also contains derivatives, 

S w (u,$, 4>) = s J tfxy/^gW™' {2uijV k Vi(t> +4 [Vi& + V^)] [V fc fe + V^)] + . . . } 

(42) 

The part (|4ip of the action gives healthy kinetic term for the (gauge fixed) vector field £ l 
and provides kinetic mixing between the fields (j) and uij. This property is precisely the 
same as in the Fierz-Pauli theory in Minkowski background [37] . The sum (Si nv + S m ) can 
be diagonalized by the shift of the field (cf. Ref. [37J), 

M 2 1 
Uy -> Uij + —r]ij(p 

As a result, the field <j> obtains healthy kinetic term 

= ~M 2 (m 2 - 2A) J d A x^g~ V<0V*0 

The overall sign here is normal (non-tachyonic) , provided that inequality (JSj) is satisfied. 

For the de Sitter or anti-de Sitter background one has Wyki = 0, so the Stiickelberg 
action is clearly second order in derivatives in that case and the fild ity describes five 
propagating modes, in agreement with Ref. [35] . 

In the general Einstein background, the Weyl tensor does not vanish, and the term (|42p 
looks dangerous, as it appears to contain four derivatives. However, integrating by parts 
and using the properties of the Weyl tensor, we write this term in the following form, 

S w (u,£, cf>) = sj d A x^W iklj {2^-VfcV^ - WWV"W</> + 4V m 0^' + 2^ j ) + ■■■} 

Hence, this term contains in fact at most two derivatives, so the number of propagating 
modes remains equal to five. The Boulwarc-Dcscr phenomenon is absent in our model, at 
least in the Einstein backgrounds. 
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In sufficiently weak background fields, when |W,yjy| <C M 2 , the mass term S m dominates 
over Sw-i so the propagating modes are not ghosts. In stronger background fields, the term 
Syy induces explicit two-derivative term in the action for the field <f>, and also extra kinetic 
mixing between this field and Uij, so there may appear ghost modes. We will comment on 
this point in section [6l 

5 Interaction between sources in Minkowski background 

Let us now consider the linearized theory in Minkowski background and study the interac- 
tion between sources. Our purpose here is twofold. First, we will confirm that all modes in 
this theory linearized about flat space-time are neither ghosts nor tachyons. Second, we will 
see that the interaction between conserved energy-momentum tensors which couple to the 
vierbein and do not directly couple to connection is mediated by both massless and massive 
spin-2 fields, so our model is indeed an infrared-modified gravity. 

Let us denote the sources coupled to the vierbein and connection by J/ 4 and S 10 ^, 
respectively, and introduce the source term in the action, 

Ssource = J d X (ih^J^ — ^AijllS'^^J (43) 

where h\ is defined by 

4 = Sj, + ^ 

We will still work with objects like J ij , S ljk , defined by J ij = J'^dj,, S ljk = S ij »5 k . Note 
that the source J lJ in general is not symmetric. 

The theory is invariant under linearized local frame rotations and linearized general 
coordinate transformations. The requirement that the source term (|43|) respects these 
gauge symmetries gives two conservation laws, 

8 j J l3 = (44) 

and 

d t S ljl = 4JM (45) 

Note that eq. (pH]) implies 

&J m = -&J M (46) 

The local symmetries also enable one to set the antisymmetric part of equal to zero. 
Thus, in what follows we use the gauge hij = hj{. 
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Let us write down the linearized field equations with the sources. We omit the subscript 
(1) in this section. The linearized gravitational equation is 



ci [Fji - ^VijF) + 9 k F ijk = J fi 



(47) 



Note that the antisymmetric part of this equation gives 



2 

which is a constraint. The propagation equation is thus the symmetric part of (|47p . namely, 

c i ( F m - + dkp mk = J(n) (49) 

The linearized torsion equations in the presence of source terms have the following form, 
c 3 (rf k d m Fi m - rf k d m F im - d i F^ k + d j F l 



+c 4 U k d m F mj - rf k d m F mi - 3 l F k i + & F 



+2c 5 (7] ik d j F - r? k &F) + Ac 6 e ijkrn d m (e ■ F) + H ijk = ^S ijk (50) 

Making use of (pQ) and ((2]) one writes the source term (|43[) in terms of hij and the components 
of the torsion. The result is 

Ssource = j d X ^2hijT ^ + — tf.[ij]S ^ + —VjS'^i — —EijkmO' S 3 ^ (51) 

where t 13 is defined by 

T a = jM) - -d m s m ^ (52) 

In view of eqs. (j4"5|) and (j4Uj) the source Tjj is conserved, as it should. 
Let us also define the following combinations of the sources, 

S = Eij kl diS ljk 

and 

<?ij = J(ij) ~ ~2 dmSm{ij) ^ 
It is now a matter of straightforward but tedious calculation to find the solution to eqs. (|47p 
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and (1501), The result is 



=hv 2 fa ~ W) ~ ^FT^ fa " W) (54) 



C\ K \ I J C\OL 

_ | gij . _|_ g C 3 

6a V 3 c\ 



S ij j + 8^ifc m a mi ) (55) 



°' = ~ 9«« 2 h 2 I 2 + 7^ £ ^ ^ + V k%k m^ (56) 

288toqC6 «r + 18a { 3a J 

i 1 r / l \ / l \ 

^ ~~ 2a A; 2 + m 2 1 \ ^ ~~ ~ j \ aik ~ ^ ik<J ) + m 2 ( k i a mj ~ kja mi/ 

36a 2 ^ i^i^jkm kj^ikm ^kk^ijrn) ^^^(.Vik^jm VjkSim ) 
- 7Tr{<%fc + \{Sikj ~ Sjki)} + 7r^k m (VikVmj ~ TljkVmi) (57) 

6a 2 3cia 

where too is the flat limit of the mass of the pseudoscalar field (see eq. (fT9|) ). 

9 a 

m. 



" " 16c 6 



and m is given by (|32p . Plugging these expressions back into the action (this amounts to 
calculating (1/2) of the source term (|5ip ) we find the action that describes the interaction of 
the sources. We write its expression omitting the terms which are ultra-local in the sources, 



<1>k ^ 1a\~,2 SS 2 + ~ 4% ~ 

' 144a k 2 + 7715 c\ k 2 J \ 2 



a 1 

aci k 2 + m 



(58) 



where bar denotes complex conjugation. The three terms here correspond to exchange by 
massive spin-0 particle, massless spin-2 particle and massive spin-2 particle, respectively. 
The latter exhibits the van Dam-Veltman-Zakharov discontinuity, as is generally the case. 
It is clear that with the restrictions on parameters summarized in , neither of the modes 
is ghost or tachyon. 

Of particular interest is the metric perturbation generated by the symmetric energy- 
momentum tensor Tjj = Tuj\ coupled to metric and not directly coupled to torsion. In this 
case we have = = Tij, = 0, drf 3 = and the expression (f54"|) becomes 

hij = 7 1 T 2 fa " fa T ) - a^^T^ fa " m 

Thus, in this case too, the interaction is mediated by both massless and massive spin-2 
fields, with relative strength being a free parameter in our model. For completeness, let us 
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write down the linearized connection in this case 




+ 



p-p— 2 Iki \ Tjk - ^VjkT^j - kj ^T ik - -7] ik 



) 



} 



It is different from the Riemannian connection corresponding to the vierbein perturbation 
(|59p (see also the first line in (|57p ). which clearly shows mixing between vierbein and torsion 
fields in our model. 

6 Conclusions 

The model discussed in this paper belongs to the class of modified gravities in the sense 
that the gravitational force is mediated by both massless and massive tensor fields. Yet 
the model successfully passes a non-trivial consistency check: in torsionless Einstein space 
backgrounds it has no pathologies in the spectrum, at least for small enough background 
curvature. Clearly, this model deserves further study. 

One issue to be understood is whether the model is consistent in more general back- 
grounds, including those with non-vanishing torsion. Another is whether the Vainshtein 
mechanism cures the van Dam-Veltman-Zakharov problem; this issue can probably be un- 
dersood in an appropriate decoupling limit, in analogy to Refs. [371 [381 El S E] • There is one 
more property of our model that may be related to the Vainshtein non-linearity. We have 
seen in section POl that at the linearized level in the Einstein backgrounds, the longitudinal 
mode 4> may become a ghost for |Wi_jfcj| > m 2 , where m is the mass of the spin-2 field. For 
spherical source, the Weyl tensor is of order |Wyjt;| ~ Rs/r 3 , where r is the distance to the 
source and Rs is the Schwarzschild radius. Hence, the danger of a ghost occurs at 



Note that is the smallest of all Vainshtein radii in the Fierz-Pauli theories [37], the 
generic value being 



By analogy to other known infrared modified gravities we expect that the longitudinal sector 
of our model goes non- linear at least at r < r%, so the Vainshtein mechanism may cure the 
potential ghost problem as well. 
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The model we studied in this paper is geometrical, and its full non-linear action is 
well defined. At the linearized level virtually every source produces perturbations of both 
metric and torsion. Thus, it would be interesting to understand whether or not black holes 
have similar property. If they do, the existence of torsionless Ricci flat solution (for zero 
cosmological constant) — Schwarzchild black hole — would mean that black holes in this 
model have torsion hair, and that the Birkhoff 's theorem is not valid. 
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